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ABSTRACT 

The relation between long gamma-ray bursts (LGRBs) and low-luminosity GRBs (ZZGRBs) is a long 
standing puzzle - on the one hand their high energy emission properties are fundamentally different, 
implying a different gamma-ray source, yet both are associated with similar supernovae of the same 
peculiar type (broad-line Ic), pointing at a similar progenitor and a similar explosion mechanism. 
Here we analyze the multi-wavelength data of the particularly well-observed SN 2006aj, associated 
with /ZGRB 060218, finding that its progenitor star is sheathed in an extended (> IOORq), low-mass 
(~ O.OIMq) envelope. This progenitor structure implies that the gamma-ray emission in this ZZGRB 
is generated by a mildly relativistic shock breakout. It also suggest a unified picture for ZZGRBs and 
LGRBs, where the key difference is the existence of an extended low-mass envelope in ZZGRBs and its 
absence in LGRBs. The same engine, which launches a relativistic jet, can drive the two explosions, 
but, while in LGRBs the ultra-relativistic jet emerges from the bare progenitor star and produces the 
observed gamma-rays, in ZZGRBs the extended envelope smothers the jet and prevents the generation 
of a large gamma-ray luminosity. Instead, the jet deposits all its energy in the envelope, driving a 
mildly relativistic shock that upon breakout produces a ZZGRB. In addition for giving a unified view 
of the two phenomena, this model provides a natural explanation to many observed properties of 
ZZGRBs. It also implies that ZZGRBs are a viable source of the observed extra-galactic diffuse neutrino 
flux and that they are promising sources for future gravitational wave detectors. 

Subject headings: gamma-ray burst: general — gamma-ray burst: individual (GRB060218, 
GRB980425, GRB031203, GRB100316D ) — supernovae: general — supernovae: 
individual (SN2006aj, SN1998bw, SN20031w, SN2010bh) - neutrinos — gravita¬ 
tional waves 


1. INTRODUCTION 

The gamma-ray emission of LGRBs and ZZGRBs show 
almost no similarities apart for being detected by the 
same instruments. LGRBs are luminous (10 50 — 10 52 
erg/s), hard (> 100 keV), highly variable and nar¬ 
rowly collimated with a typical duration of 10-100 s 
dPiranl 12004 and references therein). ZZGRBs are 
fainter by about four orders of magnitude (10 46 — 
10 48 erg/s), relatively soft (< 100 keV), not highly 
beamed and show no significant temporal variability 
over thei r entire duration, wh i ch is often longe r than 
1000 s (IKulkar ni et al I | 1998(: |Campana. et al. I 120061 
ISoderberg et al. 1120061 iKaneko et al. 1120071b 

To date there are only four well observed ZZGRBs0, 
compared to hundreds of LGRBs. However, this 
is a result of their low luminosity, which limits 
the detection to a distance of ~ 100 Mpc, com¬ 
pared to LGRBs which are seen through the en¬ 
tire Universe. In fact, ZZGRBs are much more 
common than LGRBs, and are the most abundant 
known relativistic explosions in the n earby Universe 
(jSoderberg et al. ll2006HPian et al. 1120061 ). Thus, ZZGRBs 
are of special interest, both for the understanding of 
GRBs and their connection to SNe, and as sources 
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2 These are GRBs with luminosity < 10 48 erg/s for which 
an associated SNe was observed - 11GRB/SN: 980425/1998bw, 
031203/20031w, 060218/2006aj and 100316D/2010bh. 


of high energy non-electro magneti c signals, suc h as 
gravitational-waves (e.g .. |Kotake. Takiwaki fc Harikad 

l2012tlBi rnholtz fe PirarJ2013tlAndo et al. ||2013|). neutri- 

nos (e.g.. iMurase et al. iboOftlly urase fc Iokall2013h and 
cosm ic-ravsfe.g., iBudnik et al. 1 120081 iLiu. Wang fe; Pail 

unnD- 

Based on high energy emission alone there is no reason 
to assume that LGRBs and ZZGRBs are related. More¬ 
over, theoretical considerations show that the gamma- 
rays seen in ZZGRBs cannot be produced in the same 
environ ment w h ere the gam ma-ray s in LGRBs are gen¬ 
erated (|Bromberg. Nakar fe Pira~nll2011[ ). It is therefore 
puzzling that these two apparently different GRB types 
are both associated with very s imilar peculiar SNe o f 
the rare broad-line Ic type (e.g.. iMelandri et al. 11 2 01 dll . 
These SNe show no signs of H or He in their spectra, an 
indication of highly stripped progenitors. Their ejecta 
have unusually high velocities for SNe (10,000—30,000 
krn/s), their peak luminosities indicate a relatively large 
amount of synthesized 56 Ni, and the total kinetic en- 
ergy carried by some of these SNe is unusually high 
(jWooslev k, Blooml I2006L and references therein). The 
similarity of the associated SNe suggests that ZZGRBs 
and LGRBs have similar progenitors and similar inner 
explosion mechanism. The natural question that arises 
is how similar explosions produce such different gamma- 
ray signals. 

Here we approach this puzzle by analyzing the early 
(first day) optical/UV light curve of SN 2006aj, which 
is associated with ZZGRB 060218, in order to study 
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its progenitor structure. UGRB 060218/SN 2006aj 
has the best early observational coverage out of the 
four well observed /ZGRBs and their associated SNe. 
It includes Swift continuous gamma-ray, X-ray, UV 
and optical observa tions ranging from 10 2 — 10 6 s af¬ 
ter the explosion dCampana et al. I 1200 61. many op¬ 
tical spect ra starting les s than two days aft e r the 
explosion dMirabal et alA 120061. iModiaz et al. 1 120061: 

Mazzali et al. 2006 nPian et q/*1l2006t iSollerman et al. I 

200611 a nd radio observations st arting a day after the ex¬ 
plosion (ISoderberg et al. II2006H . In fact SN 2006aj has 
probably the most detailed early optical/UV photomet¬ 
ric coverage out of the thousands SNe observed to date. 

The unique feature of the optical/UV light curve of 
SN 2006aj is that it shows two peaks. Using a recently 
develo ped method_for the analysis of double-peaked SNe 
(INakar &; Pirol I2014T) we constrain the progenitor prop¬ 
erties. These properties are then used to learn about the 
physics of 1/GRBs and on their relation to LGRBs. 

The paper structure is as follows. Section )J2] presents 
the analysis of the optical/UV light curve and the re¬ 
sulting constraints on the progenitor structure of HGRB 
060218/SN 2006aj. These constraints strongly support 
the suggestion that ZIGRBs are generated by shock break¬ 
outs (tJ31) - A unified picture for /IGRBs and LGRBs that 
naturally explain the similarities and differences between 
them is presented in jlj This picture provides a simple 
explanation to the unique velocity profile of SNe asso¬ 
ciated with 1ZGRB (!}5]). The implication of this picture 
for ZlGRBs’ neutrino and gravitational wave emission is 
discussed in IJS] 

2. THE PROGENITOR OF UGRB 060218/SN 2006AJ 

Figure Q] depicts t he optical/UV light cu rve of SN 
2006aj as taken from ICampana et al. H2006f ). It shows 
two peaks in the optical bands, at t ~ 10 hr and t ss 10 
days, where t is time since first detection of the gamma- 
rays, estimated here as the explosion time. Such double- 
peaked light curves are very rare among SNe. In typ¬ 
ical SNe the light curve is dominated by one of two 
power sources: (i) the internal energy deposited by the 
SN shock, known as “cooling envelope emission”, or (ii) 
the radioactive decay of 56 Ni. Each one of these power 
sources produces only a single peak in the optical and, 
in typical SNe, the time scales of the maximal contribu¬ 
tion to the optical light from each of the two sources are 
comparable. Therefore, observed SN light curves usu¬ 
ally contain only a single optical peak which is powered 
by the stronger power source at any given SNe. This is 
cooling envelope in explosions of red supergiants, such as 
type II-P SNe, and 56 Ni in explosions of more compact 
progenitors, such as type I and 1987-like SNe. 

Two peaks are observed in rare cases where at early 
time the emission is powered by the cooling envelope, 
which then decays sharply on a time scale comparable 
to that of the rising 56 Ni contribution. This behavior 
requires an atypical progenitor structure of a compact 
mas sive core that is engulfed by extended l ow-mass mate¬ 
rial (iHoflich. Laneer fe Dusching cr 1993; Bcr sten et al. I 
120121 INakar fc Pirol I2014H . The second peak in these 
cases is similar to the main peak of a typical 56 Ni pow¬ 
ered SNe, and thus its properties provide an estimate of 
the total ejecta mass and of the 56 Ni mass. In the case 
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Fig. 1.— Swift UVOT light curve of SN 200 6aj in five filters 
(central wavelengt h given in the l egend) from Campana et al. 

( 2006). Following [Campana et al. ( 2006) the light curve is cor¬ 
rected for Galactic reddening, E(B — V) = 0.14 (assuming a Galac¬ 
tic reddening curve with Ry = 3.1), and a host galaxy reddening 
E(B — V) = 0.20 (assuming a Small Magellanic Cloud reddening 
curve). 


of SN 2006aj the second peak is clearly powered by 56 Ni 
and it shows a total ejecta mass of ~_2 Mq, out of which 
~ 0.2 Mq are 56 Ni (jMazzali et al. Il2006f ). The only nat¬ 
ural source of the first peak in SN 2006aj is the cooling 
envelope phase of an extended mass (see appendix[A]) and 
thus its properties can provide a robust estimate of the 
radius and the mass of t he ext ended material. Here the 
results of INakar fc Pirol (|2014[ ) are used to derive these 
constraints. 

The mass of the extended material can be est i mate d 
from the time of the first peak, t p (INakar fc Pirol [2014f: 
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Vext 


0.2 c V 10 hr 
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where c is the light speed and v ex t is the velocity to which 
the extended material is accelerated by the explosion. 
Spectroscopic observations limit v ex t >0.1 c, the mea¬ 
sured photospheric velocity at t = 2.89 day dPian et al. I 
12006 ). On the high end it is most likely that v ext < 0.3 c, 
since at higher velocity the kinetic energy carried by M ext 
would be larger than the kinetic energy d eposited b y the 
explo sion in the massive core (~ 10 51 erg dMazzali et aL~\ 

l 2006 li . 

The pre-explosion radius of the extended material, 
R ex t: can be estimated by the bolometric luminosity 
at the first peak. The colors before and during the 
first peak are very blue and they are constant in time 
(jSimon. Pizzichini fc Hude5l2010f) . as expected if the ob¬ 
served bands are on the Rayleigh-Jeans part of a black- 
body spectrum during that time (i.e., with temperature 
T(t < t p ) > 50,000 K; see consistency check below). 
This implies that the total luminosity seen in UV dur¬ 
ing the peak, ~ 3 • 10 44 erg, is only a lower limit on the 
true bolometric luminosity, which may be significantly 
hig her. Since R ex t, is li near in the bolometric luminos¬ 
ity (INakar fc Pirol[20Li l. the available observations set a 
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lower limit: 

R " si0 “(-5f:)" cm (2) 

This is consistent with the lack of color evolution at 
t < t p and the model prediction that temperature 
is dropping with time, reaching at the peak T(t p ) ~ 
50,000(-R ex t/T0 13 cm) 1 / 4 K dNakar fc PirollMl . Thus, 
the optical/UV light curve of SN 2006aj indicates that its 
progenitor had a relatively compact core of several solar 
masses, surrounded by ~ 0.01 M 0 which is extended to 
a radius of a supergiant. This structure is very different 
than the typically expected structure of a fully H stripped 
progenitor, based on stellar evolution models, yet it must 
be very common in GRB progenitors given that ZZGRBs 
are more common than LGRBs. This progenitor struc¬ 
ture has several far reaching implications for the physics 
of ZZGRBs and their associated SNe, which are discussed 
in the following sections. 


3. SHOCK BREAKOUT ORIGIN FOR ZZGRBs 

The Thomson optical depth of the extended material is 
high, ~ 3,000(I? ea ; i /10 13 cm) -2 . As a result, the break¬ 
out of the shock driven by the explosion takes place at 
R ex t • Radio observations show that the leading edge o f 
the outflow is mildly relativistic ( Sodcrbcrg et al. 2006), 
implying that the breakout must be at least at a mildly 
relativistic velocity, i.e., Vbo > 0.5 c. Since rate con¬ 
siderations indicate th at the gamma-rays in Z ZGRBs are 
not strongly beamed dSoderberg et al. I [2006 1 and late 
SN spectroscopyand polarimetry show no signs of ejecta 
a-sphericity dMazzali et al. 1 120071) . the breakout is not 
expect to strongly deviate from a spherical symmetry. In 
that case the main characteristics of a mildly relativistic 
shock breakout signal, its luminosity, duration and typ¬ 
ical photon en ergy, depend only on the breakout radius 
(INakar fc Sarill2012 ): 


Uo ~ 2 • 10 


,46 
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- 1000 
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3 • 10 13 cm 
Rext 


erg s 


3 • 10 13 cm 
- 50 keV 
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This is similar to the actual gamma-ray signal of ZZGRB 


060218 where Lb 0 , 0 bs ~ 3 • 10 46 erg s t h 


1,000 


._-_i Ibo.obs/' 

s and Tbo,obs ~ 40 keV (|Kaneko et al. 1 12007 ) and it 
fits very well to a breakout radius R ex t ~ 3 • 10 13 
cm. Thus, the combination of optical/UV and radio 
observations imply that a shock breakout signal is in¬ 
evitable and that its properties are similar to the ob¬ 
served ZZGRB . As shock breakout also explains a large 
range of properties of the high energy emission from 
ZZGRBs ( e.g., smoot h profile, spectral evolution, low 
beaming: INakar fe Sarill2012l) . this result practically im¬ 
plies that the entire gamma-ray signal in ZZGRB 060218 is 
generated by a mildly relativistic shock breakout, with¬ 
out any significant contribution from a relativistic jet. 
It also lends a very strong support for the suggestion 
that all ZZ GRBs are sh ock breakouts jKulkarni_eZ_aZ 1 " 


1998 : iTan. Matzner fe McKeel [20011: ICampana, et al. 
2006: IWaxman. Meszaros fe Campanal 20071: Id 1 20071 : 


Katz Budnik fe Waxmanll2010HNakar fe Sarill2012lb 
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Fig. 2.— A schematic sketch illustrating the similarity and dif¬ 
ferences between ZZGRBs and LGRBs. Both explosions go through 
a collapse of a similar core which leads to the formation of a sim¬ 
ilar GRB engine and to a similar SN explosion. In both types 
the GRB engine launches ultra-relativistic narrowly collimated jet, 
which penetrates through the core. In LGRBs the jet is free to ex¬ 
pand as soon as it is out of the core where it produces a luminous, 
hard, narrowly collimated beam of gamma-rays which can vary in 
time on short time scales. In ZZGRB the jet emerges from the core 
into the low-mass extended material where it is choked and any 
radiation that it produces is absorbed and cannot reach to the ob¬ 
server. The jet energy is deposited in the extended material driving 
a strong shock into it. The shock is much less relativistic than the 
jet (most likely Newtonian) and it accelerates before breakout (of¬ 
ten to a mildly relativistic velocity). Upon breakout it produces 
low-luminosity soft gamma-rays which show no significant variabil¬ 
ity with time and are not narrowly beamed. 


4. A UNIFIED PICTURE FOR LGRBS AND ZZGRBs 

If all ZZGRB progenitors have a similar structure to 
that of ZZGRB 060218 then it provides a natural solution 
to the puzzle why two explosions with similar inner 
workings produce such different gamma-ray signals. 
According to the standard model for LGRBs the burst 
is powered by a central engine that launches a highly 
collimated ultra-relativistic bipolar jet. In order to pro¬ 
duce a LGRB the jet must first punch its way through 
the star and then expand freely at ultra-relativistic 
velocities to radii where generated gamma-rays can be 
seen by the observer. While the jet drills through the 
dense stellar matter its energy is dissipated and the 
engine must continue to supply power into the jet if it is 
to succeed punchi ng through th e star and produ ce the 
observed LGRB ([Z hang, Wooslev fc MacFadvenl 120031 
i Morsonv. Lazzati fe Begelmanl 120071 : iMizuta fe Alovl 
120091 iBromberg et al. 11201 111 . Thus, a necessary condi¬ 
tion for the production of a LGRB is that the engine 
working time is long enough to allow the jet to drill 
through the star. Observations indicate that a typical 
LGRB engine launches a jet at a typical isotropic 
equivalent luminosity of Li SO ~ 10 51 erg/s and a typical 
opening ang le 9j ~ 10° over a typical duration of ~ 20 s 
(jPiranl 120041) . The total energy carried by the jet, after 
correction for beaming, is ~ 10 51 erg. If the progenitor 
is a bare H stripped star of several solar masses and 
several solar radii it takes ~ 10 s for th e jet to penetrate 
throu gh the star (see appendix [HI IBromberg et al. I 
HHI), implying that the jet can successfully emerge 
from the star and that the collapse of such a progenitor 
can lead to a LGRB. 

The picture, however, is very different if there is an 
additional extended envelope surrounding the massive 

























































4 


core, similar to the one found here for ZZGRB 060218. 
Although the extended material mass is low its large ra¬ 
dius makes it very hard for a jet to penetrate. In fact 
the time that the engine must work in order for the jet to 
drill through the entire extended mass is (appendix [Bj : 


teng > 150 


( L/iso 
10 51 erg/s 


Rext \ 
3 x 10 13 cm / 


f Mext \ 

Uo- 2 M<J 


(4) 

This time is considerably longer than the typical work¬ 
ing time of a LGRB engine. Thus, a collapse of the pro¬ 
genitor of ZZGRB 060218 and the formation of a LGRB 
engine at its center will not lead to an observed LGRB. 
Instead, the launched jet, which penetrates the stellar 
core, is choked while still propagating in the extended 
material. The energy carried by the jet ( Ej et ~ 10 51 

erg) is then deposited in the extended mass accelerating 

1 1 


it to high velocity (v ex t « 0.3c 


1051 erg 


Me X t 

10~ 2 Mq 


■) 


and driving into it a strong shock. The shock accelerates 
further at the dropping density gradient near R ext and 
upon breakout produces a ZZGRB. 

Note that while the energy deposition is done by a 
narrow jet and is therefore highly aspherical, the shock 
upon breakout can be quasi-spherical. The reason is that 
the jet is choked long before it approaches R ext and the 
resulting blast wave becomes much more spherical during 
its propagation before it breaks out at R ex t ■ A schematic 
sketch of the similarities and differences between ZZGRBs 
and LGRBs according this picture is illustrated in figure 


5. THE UNCOMMON VELOCITY PROFILE OF 
SNE ASSOCIATED WITH Z/GIlBs 

This scenario resolves yet another puzzle related to 
SNe associated with ZZGRBs - why is the kinetic en¬ 
ergy in the ir fast moving ejecta is so high compared to 
other SNe dSoderberg et al. 1 120061b In typical SNe the 
explosion energy is all deposited at the center of the pro¬ 
genitor. This drives a shock that crosses first the bulk 
of the mass and then accelerates at the sharp density 
drop near the stellar edge. This acceleration dictates 
a certain relation between the kinetic energy carried by 
slow and by fast moving material, whe re E^{v) oc v~ 5 
dSakurail 119601 iMatzner fe McKee] I1999H . This relation 
is seen in regular SNe, but not in ZZGRBs where the 
fast moving ejecta carries mu ch more energy than it pre¬ 
dicts (ISoderberg et al. Il2006| h For example, in SN 2006aj 
about 2 x 10 51 e rg are carried by the bulk of the mass at 
~ 10,000 km/s (jMazzali et al. 112006t h In a regular SNe 
if a mildly relativistic (> 150,000 km/s) ejecta exist, 
the relation Ek(v) oc v~ 5 dictates that it should carry 
~ 2 x 10 45 erg. Instead, in SN 2006aj radio observa¬ 
tions indica te that the mildly relativistic material carries 
> 10 48 erg (ISoderberg et al . 1120061 iBarniol Duran et al. I 
I2014H . This observed property of SN 2006aj is naturally 
explained by the picture of ZZGRBs presented here. In 
this picture the energy in the slow moving material is 
deposited by the SN explosion mechanism at the center, 
while the observed > 10 48 erg in the fast moving material 
is deposited directly by a GRB jet, thereby decoupling 
the amount of energy carried by each of the components. 


6. NEUTRINOS AND GRAVITATIONAL WAVES 
FROM ZZGRBs 

An intriguing implication of the arising picture is the 
prospects for future detection of non-electromagnetic sig¬ 
nal from ZZGRBs. The extreme energies and velocities 
involved in LGRB engines and jets make them a poten¬ 
tial source of gravitational waves (GW), neutrinos and 
s ' high-energy cosmic rays. However LGRBs are very rare 
in the local Universe and are typically seen at a distance 
>Gpc. Here we suggest that ZZGRBs harbor the same en¬ 
gine as LGRBs, which produces similar ultra-relativistic 
jets. The propagation of an ZZGRB jet is similar to that of 
a LGRB jet within the massive core. After the jet breaks 
out of the core and into the ZZGRB extended envelope the 
envelope density is low, so the pressure in the cocoon does 
not affect the jet (see appendix [Bj . Thus, at any loca¬ 
tion that is far from the jet head the jet is unaware of the 
extended envelope. Therefore, all the physical processes 
that take place during the formation of the engine, the 
launching of the jet and the jet propagation in LGRB 
also take place in ZZGRBs up to the radius where ZZGRB 
jets are choked in their progenitors’ extended envelopes, 
namely ~ 10 12 — 10 13 cm. Thus, the same emission gen¬ 
erated by a LGRB engine and by its jet while it prop¬ 
agates up to a radius of ~ 10 12 — 10 13 cm are expected 
to be generated also by a ZZGRB . This includes photons, 
high energy particles (cosmic-rays), neutrinos and GWs. 
The extended envelope has a Thompson optical depth 
100 and therefore photons and cosmic-rays cannot es¬ 
cape through the extended envelope (the cross-section for 
pp inelastic collision at ~ TeV energies is ~ 0.1 Thom¬ 
son). However, the envelope is transparent to neutrinos 
and GWs. It is therefore worth considering the implica¬ 
tions of the model presented here for neutrino and GW 
emission from ZZGRBs, especially given that ZZGRBs are 
more common than LGRBs. 

LGRBs are expected to be bright sources of high en¬ 
ergy neutrinos (~ 10 14 — 10 16 eV). The most promis¬ 
ing production site of neutrinos is internal sho cks within 
the relativistic jet (jWaxman & Bahcallill997 ). At radii 
that are large enough > 10 11 — 10 12 cm these shocks 
are collisionless and are therefore expected to efficiently 
accelerate protons (at smaller radii the shocks are ra¬ 
diation m ediated a nd no efficient part icle acceleration is 
expected: iLevinson fc Bromberd[2008l f. At radii that are 
small enough, < 10 14 , the photon density is high enough 
to allow an efficient photo-pion production and thus a 
generation of high energy neutrinos. Current measure¬ 
ments limit the neutrino flux from LGRBs to < 

2 x I Q 10 GeV cm~ 2 s~ 4 sr ~ 4 per flavor, where Eq = 100 
TeV (lAartsen et al. Il2014f) . This upper limit is compara¬ 
ble to recent estimates of the flux expected from LGRBs 
if particles are accelerated efficien tly in internal sh ocks at 
radii of ~ 10 12 — 10 14 cm (iHummer. Baerwald fc Winter! 
I2012T) and it is about two order of magnitude lowe r than 
the measured diffuse neutrino flux (lAartsen et al. Il20l4lf . 

In the picture presented here ZZGRBs are expected to 
be a much stronger source of diffuse neutrino flux than 
LGRBfJ First ZZGRBs are more numerous. The lo- 

3 Note that in the model discussed here the neutrinos are gen- 
erated in a diff e rent environment than in the model discussed by 
IMurase &; Iokal (120131) , where they assume that ZZGRBs are gener¬ 
ated by internal shocks in a wide mildly-relativistic low-luminosity 
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cal rate of ZZGRBs , without correction for beaming, is 
~ 3 x 1CT 7 Mpc -3 yr -1 , where bea ming correct ion is 
expected to be relative ly small, < 10 dSoderberg et al. I 
120001: iPian et al. Il200fil ) . This is compared to a LGRB lo- 
cal rate, u ncorrected for beaming, of ~ 10 -9 Mpc -3 yr -1 
(|Wanderman fc Piranl I2010T) . Beaming correction in¬ 
creases the true rate of LGRBs by about two orders 
of magnitude. Thus, ZZGRBs are more frequent than 
LGRBs by about an order of magnitude. Second, The 
neutrino production of LGRBs depends on the fraction 
of high-energy protons that produce pions via interac¬ 
tion with the observed gamma-rays. This fraction, de¬ 
noted as /tt, depends strongly on the burst parameters 
and it vary between bursts. Un d er opt imal cond itions 
the estimates are < 0.1 (IWaxman fc Bahcall 1997; 

iHiimmer. Baerwald &; Winter! l2012h . In ZZGRBs how¬ 

ever, the jet is buried in an envelope that is optically 
thick to high-energy protons. Thus, energy of protons 
that in LGRBs would have been released to the host 
galaxy as cosmic-rays is converted in large part to neu¬ 
trinos via pp interactions. Thus, assuming that a large 
fraction of the LGRBs’ neutrinos are generated at radii 
smaller than ~ 10 13 cm, ZZGRBs are more efficient in 
producing ~ 100 TeV neutrinos by about two orders of 
magnitude, and possibly more (this is the product of the 
ZZGRB to LGRB rate ratio and l// w ). 

As high-energy protons are accelerate within the ultra- 
relativistic narrowly collimated jets, the neutrino signal 
is narrowly beamed as well. Since the gamma-ray emis¬ 
sion from /ZGRBs is not highly beamed, most of observed 
bursts are not expected to be accompanied by a neutrino 
signal. Thus, ZZGRBs will contribute to the diffuse flux 
but they are not suitable for a targeted point-sources 
search, similar to t he search conducted for LGRBs 
(lAartsen et al. Il2014f) . Can ZZGRBs be then the main 
source of the observed diffuse flux? lAhlers fc Halzenl 
( 2014 ) find that the sources of the diffuse neutrino flux 
produce a total energy output of ~ 10 43 erg Mpc -3 yr -1 
in ~ 100 TeV neutrinos and their volumetric rate, as¬ 
suming transient sources, must be > 10 -8 Mpc~ 3 yr -1 
(as inferred from the lack of neutrino clustering). As¬ 
suming that each ZZGRB harbor a relativistic jet with a 
typical energy of ~ 10 51 erg the total energy output in 
such jets is ~ 3 x 10 44 erg Mpc -3 yr -1 . Thus, if ~ 10% 
of this energy is converted to high-energy protons be¬ 
fore the jet is choked (i.e., at radii < 10 13 cm) then 
ZZGRBs are producing the observed diffuse flux. Assum¬ 
ing that the typical jet angle is ~ 10° the rate of ZZGRBs 
for which the neutrino beam is pointed towards Earth 
is ~ 0.5 x 10~ 8 Mpc~ 3 yr~ 4 cons istent with the limit 
derived bv lAhlers fc Halzenl (l2014f) . Thus, ZZGRBs are 
certainly viable candidates for the origin of the observed 
extra-galactic neutrino flux! 

Finally, if ZZGRBs harbor the same engine and rel¬ 
ativistic jets as L GRBs then th ey produce t he same 
GW signals (e.g. , IKotake. Takiwaki fc Harikael l2012t 
IBirnholtz fc Piranl l2013li . A difference is that while 
LGRBs are always observed close to the jet axis, the line- 
of-sight to ZZGRBs is typically away from that axis. The 
GW signal from the engine can be slightly stronger along 
the jet axis (up to a factor of 1.6 compared to an average 

jet. Our findings on the progenit or structure of HGRB 060218 
strongly disfavor the model used bv IMurase &; Iokal <|2013li 


line-of-sight observer), if its origin is q uadrup ole mass 
inhomogeneity in a rotating disk dKochanek fc Piranl 
I1993T) . Other axisymmetric GW sources in the en¬ 
gine, such as mass motions and neutrino emission, 
vanis h along t he a xis and are st rong at the equator 
(IKotake. Takiwaki fc Harikad2012h . The signal from the 
jet accelerati on is also anti-beamed an d is strongest along 
the equator (IBirnholtz fc Piranll2013l) . Thus, the off-axis 
viewing angle of typical ZZGRBs is probably an advan¬ 
tage for GW detection. The main advantage of ZZGRBs 
is their much higher rate. The volumetric rate of ZZGRBs 
is larger by about an order of magnitude than that of 
all LGRBs, including those that are unobservable since 
their gamma-ray beam points away from the Earth. If 
only LGRBs that points towards earth are considered 
then the ZZGRB rate is higher by almost three order of 
magnitude. This is important since targeted GW search 
for GRBs (e.g., following a detection of gamma-rays) is 
more sensitive than a blind search, increasing the de¬ 
tection volume by a factor of ~ 3 dKochanek fc Piranl 
I1993T ). The various predicted GW signals from the en¬ 
gine and the jet are expected to be detectable by future 
gravitational wave detectors up to a distance of ~ 100 
Mpc. The rate of LGRBs at that distance is too low to 
allow a reasonable probability for detection. However, 
about one ZZGRB take place every year within a distance 
of 100 Mpc, making it a promising GW source for future 
detectors. 

7. CONCLUSIONS 

This paper analyzes the first day optical/UV light 
curve of SN 2006aj//ZGRB 060218 finding that its pro¬ 
genitor has a compact core engulfed by an extended low- 
mass material. When the information on this structure 
is combined with the high velocities inferred from the ra¬ 
dio emitting material, it implies that the shock breakout 
form the extended material must produce a gamma-ray 
signal that is consistent with the observed ZZGRB gamma- 
ray emission. This indicates that the gamma-rays in 
2006aj/ZZGRB 060218 are generated by a mildly relativis¬ 
tic shock breakout and it strongly supports the sugges¬ 
tion that the origin of all ZZGRBs’ high energy emission 
is a shock breakout. 

These results, which are directly based on the obser¬ 
vations of SN 2006aj, naturally suggest a picture that 
unifies LGRBs and ZZGRBs, explaining how two types of 
explosions that are so different in their gamma-ray sig¬ 
nature produce very similar SNe. In this picture LGRBs 
and ZZGRBs are two manifestations of a similar core col¬ 
lapse process that leads to a similar SN explosion mecha¬ 
nism and a similar GRB central engine, where the obser¬ 
vational outcome depends only on the slight differences 
in the existence, or the lack of, a low-mass extended en¬ 
velope. This model also provides a simple explanation 
to the peculiar velocity profile seen in SNe that are as¬ 
sociated with ZZGRBs. It also implies that ZZGRBs are 
more promising sources of high energy neutrinos and 
GWs than LGRBs. ZZGRBs are viable candidates as the 
source of the observed extra-galactic diffuse neutrino flux 
and are promising GW sources for the next generation 
GW detectors. 

A final comment on the progenitor structure of SN 
2006aj. While the SN light curve constrains M ext and 
R ex t it does not strongly constrain its density profile. 
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The only requirement is that most o f th e ma ss M ext 
is concentrated near R ext rtNakar fe Pirol l2014f) . This 
is consistent with any density profile p(r) where pr 3 
increases with radius at r < R ex t and decreases at 
r > Rext- Thus, we cannot determine whether the in¬ 
ferred progenitor structure is in hydrostatic equilibrium 
or whether the extended material was thrown out to 
Rext a short time before the explosion. The latter op¬ 
tion may be more attractive given that no current stellar 
evolution model predicts a hydrostatic structure simi¬ 
lar to that SN 2006aj for a fully H stripped star, while 
recently there are several lines of evidence that massive 
stars go through a strongly enha nced mass-loss episode a 
short time before they explode (IQfek et al. 112013112014 


iGal-Yam et al. 112014 iSvirski &; NakaH 120141) . Yet, an¬ 
other intriguing speculation is that the progenitor is af¬ 
fected by a binary, or maybe even by a binary merger 
(e.g .. along similar lines to those suggested bv iChevalierl 
12012f) . that put the extended material at place a short 
time before the explosion. 
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APPENDIX 

THE POWER SOURCE OF THE FIRST OPTICAL PEAK 

The analysis of the progenitor structure of SN 2006aj relies on the identification of the first optical/UV peak as a 
cooling envelope emission. Here this identification is justified by considering known and speculated emission power 
sources in SNe and GRBs. The conclusio n is that while cooling envelope emission provides a natural explanation for 
the first peak (as discussed in length in lNakar fe Pirol[2014f) . all other sources are either ruled out or are highly unlikely. 

Cooling envelope: The energy source of cooling envelope emission is the shock that crosses the star and any surround¬ 
ing mass, if it exists, in regions where the diffusion time is longer than the expansion dynamical time. In these regions 
the internal energy deposited by the shock is trapped by the gas and it cools adiabatically during the gas expansion, 
hence the term “cooling envelope”. As the outflow expands its optical depth drop and so does the radiation diffusion 
time to the observer, while the expansion time grows. In regions w here the two time scales become comparable the 
radiation escapes to the observer. As discussed in iNakar fe Pirol (|2014f) this source of emission provides a natural 
explanation for the first optical peak observed on a time scale of ~day in case of a double peaked SNe. Important 
supporting evidence for that in the case of SN 2006aj is the very blue color of the first peak and the fact that the 
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colors do not vary with time, which indicates that the band that we observe are most likely on the Rayleigh-Jeans 
part of a blackbody spectrum. This is expected for cooling envelope emission, where the optical depth at the source 
is high and the radiation has enough ti me to a chieve thermal equilibrium before it escapes, even if the outflow is fast, 
after the ejecta expanded considerably rtNakar fo Saril 12010 ). 

Radioactive decay of 56 Ni: This power source of energy, which dominates many SNe light curves, deposits energy 
at the known decay rate of 56 Ni, first to 56 Co and then to 56 Fe. The observed luminosity form 56 Ni is the energy 
deposited by radioactive decay in “exposed” regions, from where radiation can escape over a dynamical time. Since 
the total amount of mass in exposed regions depends on time and velo city, the ev olution of luminosity generated by 
56 Ni for a given outflow is set by the fraction of 56 Ni in that region dPiro fe Nakarll2013f) . This sets, at any given 
time, a maximal luminosity that 56 Ni can produce, which is the luminosity of an outflow that is composed purly of 
56 Ni. The exposed mass at the first peak is ~ 0.01(u/0.2c)Mq implying that the maximal contribution of 56 Ni to 
the luminosity at this time is ~ 10 42 (i>/0.2c) erg/s. This rules out 56 Ni as the source of the first peak which show a 
luminosity > 3 x 10 44 erg/s. 

Interaction (afterglow): Another power source seen in some SNe is a continuous interaction with the circum-stellar 
medium. In GRBs such interaction is the source of the afterglow. The difference between continuous interaction and 
cooling envelope emission is that in the former the shock takes place in a region with optical depth that allows for 
radiation to escape immediately over a dynamical time scale. Thus, if the first peak is generated by interaction then 
its luminosity is limited by the instantaneous strength of the interaction. Namely, the explosion ejecta must drive a 
strong shock into the circum-stellar medium at least up to t « 10 hr, at which point either the interaction dies (e.g., 
due to a sharp drop in the circum-stellar density) or the shock becomes radiatively inefficient. The radio emission at 
t = 1.89 day is presumably genera ted by s uch interaction and it shows that the interaction shock is propagating at a 
velocity close to the speed of light (ISoderberg et aZ.ll2006f) . The interaction at this point is much too weak to account 
for the optical emission at this epoch, but assuming that at t ss 10 hr the interaction have been much stronger, could it 
then be the source of the optical/UV? Considering all the outcomes of the entire allowed phase space for interaction is 
beyond the scope of this paper, however several general considerations show that it is highly unlikely that interaction 
can produce the observed first optical/UV peak for two reasons - it predicts an optical/UV spectrum that is too red 
and an X-ray flux that is too bright compared to the observations. 

As the interaction shock is mildly relativistic its radius at the first peak is r ~ 10 15 cm. The circum-stellar medium 
must be optically thin for Thomson scattering at this radius, otherwise the mildly relativistic shock breakout would 
have been taken place at ~ 10 15 cm, resulting in a much brighter and longer signal in gamma-rays then observed 
(equation [3]). Emission from optically thin mildly relativistic shocks are seen in late stages of GRBs and in some 
SNe. In these cases the shock is collisionless and it generates strong magnetic fields and accelerates electrons to a 
power-law distribution. As a result the dominant emission is synchrotron and the spectrum above the self absorption 
frequency (which is typically in the radio or mm bands) is a power-law that is spread over many orders of magnitude 
in frequency with a specific flux F v oc v a with a ~ — 1. This is very different than the observed UV spectrum where 
a ~ 2, which requires the self-absorption frequency to be > 10 15 Hz. However, even the highest possible circum-stellar 
density that is optically thin for Thomson scattering at ~ 10 15 cm does not bring the synchrotron (or free-free) self¬ 
absorption frequency of a mildly relativistic shock into the UV. In addition, the synchrotron power-law spectrum also 
predicts an X-ray luminosity that is comparable or larger than the UV luminosity, regardless of the location of the 
self-absorption frequency. In reality at the time of the first peak the X-ray luminosity is fainter than the UV by two 
orders of magnitude. 

Continuous central engine activity: The last power source that is often considered in GRBs and sometimes also 
in SNe is a continuous energy supply by a central engine, possibly an accreting black-hole or a long lived magnetar. 
The existence of such sources in SNe is still rather hypothetical, while in GRBs there is stronger evidence that the 
central engine can be active also on time scales of hour or days. Nevertheless, the optical emission is highly unlikely 
to be powered this way. The reason is that the bulk of the SN ejecta mass, ~ 2M 0 , lies between the center of the 
explosion and the observer. The photons observed in the first optical peak must be generated at larger radius than 
that of the bulk of the ejecta. If the energy from the central engine is deposited first in the ~ 2M 0 ejecta it is 
converted to kinetic and thermal energy of the ejecta before radiated away after the ejecta optical depth drops, over 
time scale of weeks. Thus, similarly to LGRBs, the energy generated by the central engine must “penetrate” through 
the bulk of the mass before being dissipated to optical/UV photons. Again, like in LGRBs, this may be done if the 
engine continuously launches relativistic jets. However, based on GRB observations, the expected optical/UV/X-ray 
emission from relativistic jets suffers from the same problems of interaction emission. It usually show a power-law 
spectrum with a <C 2, which does not fit the observed optical/UV spectrum and the faint X-ray emission. More 
importantly, a relativistic jet must open a cavity in the SN ejecta inducing strong spherical asymmetry in the slow 
movin g materi al, whic h is ruled out by the lack of polarization and by the spectral line profiles seen in the SN nebular 
phase ( Mazzali et al. Il2007h . 

JET PROPAGATION IN THE CORE AND IN THE EXTENDED MATERIAL 

The general phy sics of a relativistic hydrodynamical jet that propagates in a surrounding medium is described at 
length in lBromberg et ai] ( 2011 ). Here we provide a brief outline of this system, focusing on the time that the engine 
must work for the jet to penetrate through the core and through the extended material. A jet that propagates in 
surrounding media forms a forward-reverse shock structure at its head. The high pressure plasma in the jet head spills 















sideways continuously as the jet propagates forming a hot cocoon that engulfs the jet. This cocoon may or may not 
collimate the jet, depending on the jet and the external medium properties. Since energy is leaving the jet head into 
the cocoon continuously, the head propagation depends on a continuous supply of energy, which is injected into the 
head by the jet via the reverse shock. Thus, in order for the jet to propagate a given distance the engine that launches 
the jet must work for a duration that is long enough so a fresh jet material will continue to cross the reverse shock 
during the entire propagation. Thus, if the head propagates up to a distance r at a velocity f3hC the jet working time 
must be: 


teng{r) 


^ (1 -« 


(Bl) 


where the term 1 — fih includes the relative velocity between the relativistic jet and the head. This term is ~ 1 for a 
Newtonian head, implying that the engine working time is simply the jet propagation time. However, if the head is 
relativistic then by the time that the engine stops working the head is at a distance ~ ct en g from the center and the 
last jet element that was launched by the engine will catch up with the head only after t eng /{ 1 — /3h)- During that 
time the jet will continue to drive the head forward. Thus, the engine working time needed for a relativistic head to 
propagate a distance r is much shorter than r/c. 

The evolution of the jet is determined by finding the properties of the various components (e.g., head, cocoon, etc.) 
of that system. The propagation velocity of the head is set by the balance of the jet luminosity per unit area into the 
head and the ram pressure of the ambient medium. It is therefore useful to define a dimensionless parame ter w hich 
is t he ratio between the energy density of the jet and the rest-mass energy density of the ambient medium (IMatznerl 
12001 
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where Lj is the total jet luminosity, E j is the jet cross-section at the head and p is the ambient medium density at 
the head location. The propagation velocity of the jet head is: 
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Thus, the head is relativistic when L 1 and Newtonian for L < Cl. The collimation of the jet depends also on the 
half opening angle at which the jet is launched, 6 0 , where for L < d Q 4 ^ 3 the jet is collimated by the cocoon pressure. 
The jet collimation affects the value of E j and thus also the value of L. For a given set of parameters IBromberg et al. I 
( 2011 ) obtain: 
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(B4) 


where t is the time since the jet launching started and p is the external density near the jet head location. Equations 
IB3I and [B4I together can be solved to find the jet location at time t. 

The isotropic equivalent luminosity of a typical GRB jet is L iso ~ 10 51 erg/s and its opening angle is 9 0 ~ 10°. The 
beaming corrected jet luminosity is then Lj = Li SO 0Q/2. In a massive (M core > Mq) and compact ( R CO re ~ 10 11 cm) 
core the density is high and L < 1, resulting in a Newtonian (or at most a mildly relativistic) collimated jet. Thus, 
the engine working time must be comparable to the time needed for the jet to cross the core: 


7 ( L lso \ ~ 1/3 ( e 0 V /3 f R core V /3 ( M core V /3 
eng,core y IQ 51 erg s _1 / \ 10 deg/ \ 10 11 cm/ \10 Mq/ 


(B5) 


Here we used the approximation (3h « L 1 / 2 which is appropriate for Newtonian heads. 

The density of the extended material is much lower than in the core. As a result, for a typical GRB jet L > 9 Q 4 ^ 3 , 
resulting in an uncollimated relativistic jet. The minimal engine working time for a successful jet penetration is shorter 
than the extended material light crossing time, but it is still much longer than the time it takes the jet to cross the 
core: 
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where we used the approximation for a relativistic head 7 h 


L 1 / 4 /\/2, where 7 h is the head Lorentz factor. 





















